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Abstract 



Based on the ACV approach to transplanckian energies, the reduced- action 
model for the gravitational 5-matrix predicts a critical impact parameter b c ~ R = 



CI 



2Gy / s such that 5-matrix unitarity is satisfied in the perturbative region b > b, 
while it is exponentially suppressed with respect to s in the region b < b c that we 
think corresponds to gravitational collapse. Here we definitely confirm this state- 
ment by a detailed analysis of both the critical region b ~ b c and of further possible 
contributions due to quantum transitions for b < b c . We point out, however, that 
the subcritical unitarity suppression is basically due to the boundary condition 
which insures that the solutions of the model be ultraviolet-safe. As an alternative, 
relaxing such condition leads to solutions which carry short-distance singularities 
presumably regularized by the string. We suggest that through such solutions — 
depending on the detailed dynamics at the string scale — the lost probability may 
' be recovered. 
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1 Introduction 

X. 

Interest in the gravitational S-matrix at transplanckian energies [HH] has revived in the 
past few-years [SHI], when explicit solutions of the so-called reduced-action model [I] have 
been found [5]. The model is a much simplified version of the ACV eikonal approach [TJ[3] 
to transplanckian scattering in string-gravity, and is valid in the regime in which the 
gravitational radius R = 2Gy/s is much larger than the string length X s = y/a'h, so that 
string-effects are supposed to be small. 

The reduced- action model (sec. was derived by justifying the eikonal form of the 
S'-matrix at impact parameter b on the basis of string dynamics and by then calculating 
the eikonal itself (of order ~ ^ 3> 1) in the form of a 2-dimensional action, whose power 
series in ^ corresponds to an infinite sum of proper irreducible diagrams (the "multi- 
H" diagrams [2111]), evaluated in the high-energy limit. The model admits a quantum 
generalization [6] of the S'-matrix in the form of a path-integral — with definite boundary 
conditions — of the reduced-action exponential itself. 



1 



The main feature of the model and of its boundary conditions is the existence of 
a critical impact parameter b c ~ R such that, for b > b c the S-matrix matches the 
perturbative series and is unitary, while for b < b c the field solutions are complex-valued 
and the elastic S-matrix is suppressed exponentially. The suppression exponent is of order 
^ ~ §r (Ap being the Planck length) or, if we wish, of the same order as the entropy of 
a black-hole of radius R. From various arguments we believe that in the region in which 
b < b c (that is, b is smaller than the gravitational radius), a classical gravitational collapse 
is taking place. 

The model, in its simplest axisymmetric form, is formulated in terms of only one 
effective field p(r 2 ) depending on the transverse radius squared r 2 = x 2 and is defined by 
p(r 2 ) = r 2 (l — (27ri?) 2 ^2 0), where h = V 2 determines the transverse gravitational field 
and the corresponding metric, which is of shock-wave type. A key role in the derivation of 
the above features is played by the boundary condition p(0) = which avoids a possible 
singularity of h at r 2 = and is the main cause of the suppression of the S-matrix for 
b < b c . In fact, the complex solutions of the semiclassical approach for b < b c were 
interpreted at quantum level [6] as due to a tunnel effect, required in order to reach p = 
at r 2 = 0, across a potential barrier occurring in the lagrangian. 

At this point the question was (and is): do we reach S- matrix unitarity for b < b c by 
summing over inelastic processes? do we recover full information [8] from the scattering 
experiment in the collapse region? It was already found in [7J, by semiclassical methods, 
that this is not the case, and that the unitarity defect persists when all inelastic channels 
are included, although in a way dependent on the rapidity phase space parameter Y. 
This result is puzzling because one would like to know where does the probability go if 
the model is complete or — if it is not — how to complete it. 

The purpose of the present paper is both to look at possible flaws in the result just 
quoted and to suggest a tentative answer to the ensuing question. We exclude flaws in 
two ways. In sec. [3] we perform a detailed analysis of the solutions of the semiclassical 
unitarity equations, we choose the stable ones and we investigate the unitarity behaviour 
by a perturbative method around the critical point b ~ b c and by numerical methods 
elsewhere. There are no surprises: the results of [7] are fully checked, we only gain some 
better understanding of the unitarity defect around the critical point. 

As a second attempt, we look for a quantum treatment of inelastic S'-matrix ele- 
ments (sec. H]). Since the quantum version of the reduced action model features a quan- 
tum mechanical hamiltonian with a Coulomb potential in p-space, we eventually evaluate 
quantum transitions from the basic tunneling wave function to other states of the sys- 
tem. Of particular interest are the bound states in the strong-coupling region p < 
where the potential is attractive, because they could correspond to collapsed matter. Un- 
fortunately, all relevant matrix elements carry the same exponential suppression as the 
tunneling amplitude itself, and the unitarity defect survives. 

As a final point, in sec. |5]we test the boundary conditions of the model, by letting p(0) 
fluctuate away from zero, with the weight assigned to it by the reduced action itself. We 
find that, while the elastic S'-matrix element is stable — that is, dominated by very small 
p(0) - - large inelastic contributions may come from the short-distance region, where 
however the model is inadequate and string effects are expected to play an important 
role. We argue on this basis that in the direction of such ultraviolet-sensitive solutions 
the model is incomplete and that, by completing it with the proper string dynamics we 
may discover where the probability goes. 
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2 The reduced-action model for gravitational ^-matrix 



We provide here a brief account of the model being considered — as develope in refs. [SIE] 
- for both completeness sake and in order to emphasize some points which are useful in 
the following. 



2.1 The semiclassical ACV results 

The simplified ACV approach [5] to transplanckian scattering in the regime R = 2Gy / i ^> 
X s is based on two main points. Firstly, the gravitational field g^ u = rj^ u + h^ u associated 
to the high-energy scattering of light particles, reduces to a shock-wave configuration of 
the form 

h-\ x+=0 = (2rrR)a(x)6(x-) , h ++ \^ =Q = (2ir R)a(x) 6 (x + ) (la) 
hij = (7tR) 2 Q(x + x-) (sa - h(x) , (lb) 

where a, a are longitudinal profile functions, and h(x) = V 2 is a scalar field describ- 
ing one emitted-graviton polarization (the other, related to soft graviton radiation, is 
negligible in an axisymmetric configuration). 

Secondly, the high-energy dynamics itself is summarized in the /i-field emission-current 
T-L(x) generated by the external sources coupled to the longitudinal fields a and a. Such 
a vertex has been calculated long ago [HIHD] and takes the form 

- V 2 U = V 2 aV 2 a - V^-aViV^-a , (2) 

which is the basis for the gravitational effective action [TTHT3"] from which the shock-wave 
solution (pQ) emerges |4j. It is directly coupled to the field h and, indirectly, to the external 
sources s and s in the reduced 2-dimensional action 



A 



2nGs 



d 2 * 



1_ ^_ (ttR) 
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as + as- -VaVa + (-(V 2 0) 2 - 2V0 ■ VH) 



(3) 



which is the basic ingredient of the ACV simplified treatment. Note that here the gravi- 
tational radius R plays the role of (dimensionful) coupling constant and that — because 
of the higher derivatives of <ft involved — non-renormalizable UV divergences may occur 
in general. 

The equations of motion (EOM) induced by (jHJ) provide, with proper boundary condi- 
tions, some well-defined effective metric fields a and h which are, hopefully, UV-safe. The 
"on-shell" action A(b,s), evaluated on such fields, provides directly the elastic S- matrix 

S el = exp (~A(b, s) j . (4) 



Then, it can be shown jUE] that the reduced-action above (where now R/b plays the role 
of effective coupling constant) resums the so-called multi-H diagrams (fig. [T]), contributing 
a series of corrections ~ (R 2 /b 2 ) n to the leading eikonal, as well as their resummation for 
R/b = Oil). 
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Figure 1: Diagrammatic series of H and multi-H diagrams. 



Furthermore, the S'-matrix fl4]) can be extended to inelastic processes on the basis of 
the same emitted-graviton field h(x). In the eikonal formulation the inelastic S'-matrix is 
approximate!}!!] described by the coherent state operator 



S = exp (^A{b, s)^ exp \A2irR-\fai J d 2 x h(x)Q(x] 
d 2 k dk 3 



a 



Gs 



[A(x),A\x')] = Y5(x-x 



[a(k, k 3 )e ik x + h.c] = A(x) + A\x) 



(5) 
(6) 



where the operator Q(x) incorporates both emission and absorption of the /i-fields and Y 
parameterizes the rapidity phase space which is effectively allowed for the production of 
light particles (e.g. gravitons). 

For a given value of the "gravitational coupling" a = Gs/h the parameter Y is possi- 
bly large for large impact parameters b 3> V GH, because the effective transverse mass of 
the light particles is expected to be of order h/b, i.e., much smaller than the Planck mass, 
thus yielding roughly Y ~ log(sb 2 /h 2 ) ^> 1. On the other hand, we should notice that 
dynamical arguments based on energy conservation [T4] and on absorptive corrections of 
eikonal type, consistent with the AGK cutting rules [TJ] , tend to suppress the fragmenta- 
tion region in a 6-dependent way, so as to constrain Y to be O (1) for impact parameters 
in the classical collapse region b = 0(R). However, such arguments do not take into 
account possible dynamical correlations coming from multi-H diagrams, as mentioned in 
footnote [TJ It is fair to state that a full dynamical understanding of the Y parameter is 
not available yet, and for this reason we shall discuss what happens for any values of Y. 

In the case of axisymmetric solutions, where a = a(r 2 ), a = a(r 2 ), <fi = <p(r 2 ) it is 
straightforward to see, by using eq. (T5J), that l-i{r 2 ) = (d/dr 2 )'H(r 2 ) = — 2aa becomes 
proportional to the a, a kinetic term. Therefore, the action ([3]) can be rewritten in the 
more compact one-dimensional form 



A 



2n 2 Gs 



dr< 



a(r 2 )s(r 2 ) + a(r 2 )s(r 2 ) — 2paa — 



{2itR) 



■ - d(1 (7\ 



where we have introduced the auxiliary field p(r 2 ) 

p = r 2 (1 - (2tt j R) 2 0) , h = 4(r 2 0> 



P) 



1 The coherent state describes uncorrelated emission (apart from momentum conservation [2]). How- 
ever, the eikonal approach based on eq. (J3]) also predicts [S] correlated particle emission, which is sup- 
pressed by a power of (Gs/h)Y relative to the uncorrelated one, and is not considered here. 
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which incorporates the ^-dependent interaction. The external sources s(r 2 ), s(r 2 ) are 
assumed to be axisymmetric also, and are able to approximately describe the particle- 
particle case by setting ns{r 2 ) = S(r 2 ), r Ks{r 2 ) = 5(r 2 — b 2 ), where the azimuthal averaging 
procedure of ACV is assumedH 

The equations of motion, specialized to the case of particles at impact parameter b 
have the form 



2vrp 



9 



2p 2 



9(r 2 



2vrp 
(r > b) 



9(r 2 - b 2 ) 



(9) 
(10) 



and show a "Coulomb" potential in p-space, which is repulsive for p > 0, acts for r > b 
and plays an important role in the tunneling phenomenon. By replacing the EOM ([9]) 
into eq. (J^J), the reduced action can be expressed in terms of the p field only, and takes 
the simple form 



A 



-Gs / dr 



R 2 



(1 - pf - -0(r 2 - b 2 ) 



dr 2 L(p,p,r 2 ) 



111 



which is the one we shall consider at quantum level in the following. 

The effective metric ([1]) generated by the axisymmetric fields p, a and a is calculated [5] 
on the basis of the complete form of the shock-wave ([T]) and is given by 

ds 2 = -dx+daT [I - ie(x + oT)(l - p)] 

+ (dx + ) 2 5(x + ) [2irRa(r 2 ) - 1(1 - p)\x'\] 
+ (dx-) 2 5(x-) [2nRa(r 2 ) - 1(1 - p)\x + \] 



+ dr 2 



1 + 2(nR) 2 Q(x + x- 



d6 2 r 2 



l + 2(7ri?) 2 (0 + r 2 



(12) 



This metric is dynamically generated and may be regular or singular at short distances, 
depending on the behaviour of the field solutions themselves. 



2.2 Boundary conditions and critical impact parameter 

Two boundary conditions are necessary to solve the equation of motion (Tl0|) : first of all we 
set p(oo) = 1 in order to have a gravitational field h ~ 1 — p vanishing at large distances. 
The second boundary condition is p(0) = and it is necessary in order to obtain UV-safe 
solutions at short distances. Indeed, by the definition of p — which embodies an r 2 factor 
in eq. (JBJ) — , the fields (ft and h are singular if p(0) ^ 0. More precisely <j) ~ — p(0)/r 2 
is quadratically divergent, and this implies that the outgoing flux of V0 at the origin 
is ~ — p(0) and therefore h = V 2 ~ — p(0)5(r 2 ) is singular too. Thus, solutions with 
p(0) possess UV singularities, and this implies in turn an UV divergent action and a 
singular metric. In particular, the singularity of is present in the metric coefficient h rr , 
which changes sign also, while a S(r 2 ) singularity of difficult interpretation occurs in the 

2 The most direct interpretation of this configuration is the scattering of a particle off a ring-shaped 
null matter distribution, which is approximately equivalent to the particle-particle case by azimuthal 
averaging [5]. 



5 



rest of the metric (fl2|) . For these reasons we set p(0) = so as to obtain a finite action 
and a regular metric. 

The solution of the equation of motion satisfying the above boundary conditions is 

t b r 2 (r 2 < b 2 ) 

R 2 cosh 2 x(r 2 ) (r 2 >6 2 ) (13) 
r 2 — b 2 = R 2 (x + si n h X cos h X — Xfo — srn h Xfe cos h X&) 

with tanhxh = h = p{b 2 ) and Xb = x(^ 2 )- The parameter £ b is determined by requiring 
regular matching of the solution at r 2 = b 2 , that is 

t b {l-t 2 b ) = ^ (14) 

This equation acquires the meaning of criticality equation. In fact it can be solvecH only 
if the impact parameter b exceeds a critical value b c given by b 2 = ^^-R 2 . For b > b c , one 
of the two real solutions matches the perturbative result [3JI5] 

^(M = G S (log^ + |^ + --.) , (15) 

where L is an infrared cutoff parameterizing the infinite "Newtonian" phase. If instead b < 
b c , the criticality equation has two complex conjugate solutions, providing an imaginary 
part of p in eq. ( FT3|) that we can't interpret at purely classical level. Anyway we can try 
to use these complex solutions and substitute them in the on shell action, which is easily 
calculated to be 



A = Gs! [ 2[ X (L 

where x{L 2 ) — \og(L/R). 



1 -Xb]- 1 —^} (16) 



For b < b c this expression becomes complex-valued, providing a non unitary S = e 1 " 4 . 
It can be shown by stability arguments [5] that the physical solution of (JHj), when b < b c , 
is the one with negative imaginary part of t&, giving rise to a suppression of the S-matrix in 
the elastic channel. For example, for b = 0, we obtain Xb=o — — ivr/2 and the exponential 
suppression 

S ~ e~ wGs . (17) 

The suppression exponent Gs ~ R 2 /\ 2 P is of the same order as the entropy of a black- 
hole of radius R, but occurs here because of the complex-valued p{r 2 ) when the impact 
parameter is smaller than the critical value. 

We conclude this subsection by noting that real-valued solutions to eq. ffTOj) and sat- 
isfying p(oo) = 1 are always of the form (Tl3~J) and have necessarily p(0) > whenever 
b < b c . The minimum (and closest to 0) value of p(0) is found for a particular value of 
the initial slope p(0) = t m determined by 



— = cosh ( Xm ) sinh(x m ) = ^ _ f2 y2 



m) 



The corresponding solution provides, so to say, the real solution with smallest "distance" 
to the complex solution having p(0) = 0. 



3 We consider only solutions with tb > 0, otherwise the condition p(oo) — 1 would not be satisfied. 
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2.3 The quantized ^-matrix and tunnel effect 



The suppression of the S'-matrix for b < b c can be interpreted by generalizing the scat- 
tering matrix to a quantum level. In the CC proposal [6] this is achieved by summing 
the reduced-action exponential over every path p(r 2 ) satisfying the boundary conditions 
p(0) = and p(oo) = 1 



S(b,s) 



p(oo)=l 



p(0)=0 



[Vp(r 2 )} e _i 1 L ^ r2 ) Ar2 e 2J ^ N-AW") ^ 



(19) 



where we have included the coherent state operator Q of eq. ([6]) in order to describe 
inelastic processes. The S'-matrix in the elastic channel is obtained by evaluating the 
vacuum expectation value and we obtain 



Li! 



rp(oo)= 


Jp(0) 


=0 


1 




AG 


(1 



[Dp] e 



-i/L a (r 2 )dr 2 



iy)(i - p) 2 - y6(r 2 - b 2 ) 



2Y 



y = 



71 



(20) 
(21) 



The parameter y is related to the rapidity phase space allowed for the emitted gravitons, 
and provides a related absorption of the elastic channel. In this section, we will examine 
only the case y = and L y — > L, neglecting such absorptive effects. In this way, the 
model is less complicated, nevertheless it explains the suppression (fTTj) . The general case 
y 7^ has been discussed in [7]. 

The problem of calculating S e \, by use of the Trotter formula, turns out to be equiva- 
lent to quantize a hamiltonian and to evaluate a matrix element of the evolution operator 
U(0, oo). In other words, the classical dynamics of the field p(r 2 ) governed by the la- 
grangian ( II ip . is promoted to a one- dimensional quantum system where p plays the role 
of "position variable" and r 2 = r represents the "time" evolution variable. The classical 
hamiltonian is given by the Legendre transform of L and we quantize it by imposing canon- 



ical commutation relations: introducing the conjugate momentum II = |p 
we have 



1 

AG 



R 2 



P 



H = Up — L 

[p,n] = ih 

The Hamiltonian, according to this quantization, is 

'0(r 



P 

P ~ l 2a dp' 



0(r - b 2 



r = r 



Mp- 1 )' 

(22) 
(23) 



H 
1 



Aa dp A 



+ a 



P 



1 

R 1 



H 9(r - 

~z~ + a 

n p 



In this way, the path integral (|2"U|) is equivalent to the matrix element 

S(b, s) = (p = 0\U(0, +oo)|n = 0) 



(24) 



(25) 



where the initial and final states express the boundary conditions for p (we recall the 
relation | H = 0} = \p — 1}). We note that the Hamiltonian is characterized by a Coulomb 
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potential that is attractive for p < and repulsive in the region p > 0: there is an infinite 
Coulomb barrier separating the boundary condition p = from the (perturbative) region 
of large p > and by means of this feature we'll be able to interpret the suppression of 
the S- matrix as a tunnel effect. 

The quantum model can be solved for any 6, but is particularly simple for 6 = 0, 
where it explains the value S ~ q-^Gs obtained at semiclassical level. Indeed, if the 
impact parameter vanishes, the expression (|25p gives 

S(b = 0,s) = (p = 0|W c (0,+oo)|n = 0) (26) 

where U c is the evolution operator related to the time independent Coulomb hamiltonian 
H c = H + Gs/p. We calculate the above matrix element by introducing the Coulomb 
wave function 

V> c (p) = (p|W c (0,oo)|II = 0) (27) 

so that S(b — 0, s) — ?/> c (0). Therefore, the 6 = quantum problem is very similar to the 
calculation of the Coulomb wave function at the origin for nuclear processes. It can be 
shown [7] that, for every 6, the 5-matrix is given by the formula 
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S(b, °) = ^J d P e- iQ( ^ +&2) MP) (28) 

which was treated in full details in refs. [6JI7]. 

In order to derive the function ip c (p), we note that the state | H = 0) is an eigenstate 
of the free hamiltonian H with zero energy, so that \i/j c ) becomes an eigenstate of the 
full hamiltonian with null eigenvalue. Therefore, the wave function can be determined by 
solving the stationary Schrodinger equation with zero energy (from now on we use R — 1 
as unit length) 

1 d 2 ^ c , fl 



H c ip c = h 



+ a - - 1 V. 



. (29) 



4a dp 2 \p 

The form of ip c is specified, including its boundary conditions, by the Lippman-Schwinger 
equation 

ipc = e 2[ap + aGo(0) pv MP) (30) 

with G(E) = [E — Ho — ie]^ 1 . @ For p > 0, where the potential is repulsive, the Coulomb 
function contains incident and reflected waves, while it has only a transmitted wave in 
the region p < of attractive potential. This explains why the wave function ip c (p) 
is suppressed in the origin by the tunnel effect through the Coulomb barrier and the 
tunneling amplitude gives the order of magnitude of the S-matrix for 6 = 0. 

Indeed, by solving eq. (1291 with the above boundary conditions, the wave function ip c 
turns out to be [7] 



2\ia 



UiaL . 

ipc = T7 r ze 2 

cosh(7ra) 



i7T©(i,2) 

U(l + ia, 2, z) + . v .7 1 + ia, 2, z) 
1 (ia) 



Aiap) (31) 



J The — ie prescription is related to the r 2 -antiordering of W c (0,oo) 
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where L is an infrared cutoff that regularizes the Coulomb singularity at large distances 
while U and F are the irregular and regular confluent hypergeometric functions respec- 
tively. The value at the origin of the wave function is 



tl> c (0) = S(0,s) 



(AaL 



T(l + \a) cosh(7ro;" 



O 



(32) 



In this way we obtain an interpretation for the suppression of the elastic amplitude for 
b < b c , as well as quantum corrections to it. 



3 Semiclassical unitarity defect below the critical point 

The reduced-action model described in the previous section shows the existence of a 
critical impact parameter b c ~ R(s) such that, for b < b c the elastic scattering amplitude 
is exponentially suppressed with s both at semiclassical and quantum level, even without 
inelastic processes (y = 0). For b > b c the small elastic absorption which appears in the 
quantum model at y > is expected to be compensated by emission processes. We do 
not know about b < b c , but it might be possible that the elastic suppression due to the 
tunneling is still compensated, because of the quantum gravity dynamics. 

In any case, the model can be unitary only if we consider the total transition proba- 
bility, summed over all emission processes. For instance, at b — 0, the suppression factor 

|(0|S|0)| 2 ~e- 2 ™ (33) 

could be compensated, thus recovering unitarity, only in the sum of all probabilities 

J2\(n\S\0}\ 2 =(0\S^S\0) = l, (34) 

n 

where \n) generically indicates "the" state with n gravitons which can be emitted thanks 
to the coherent state operator in eqs. fl6|T9|) . 



3.1 Inclusive action and equations of motion 



In order to study the unitarity properties of our model, we try to evaluate the vacuum 
expectation value of S^S. From the definition ( TT9|) . it is clear that the result will be a 
double path-integral in the fields p(r) [from S] and p{r) [from S*]. The matrix element 
between the vacuum states is computed after normal ordering of the operators A, A< in 
the coherent-state operators by means of the Baker-Campbell-Hausdorff formula 



dr 



We obtain the expression 



(0|Ste|0)= [Vp][Dp}e 



(35) 



(36) 



where the usual boundary conditions p(0) = p(0) = 0, p(oo) = p(oo) = 1 are understood 
and we introduced the inclusive action 

b 2 ) ^ &(T-b 2 ) 



A u 



-a 



py 



e( 



T 



P 



(1-P) 2 + 



~\2 



P 



iy(p - p) 



dr . (37) 
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Here the dependent term shows the effect due to summing over ail intermediate states. 

In the spirit of the semiclassical approximation, a good estimate of the value of the 
functional integral for large a, i.e., for transplanckian energies, is obtained by evaluating 
the inclusive action on the field configuration that maximizes iA u . The inclusive action 
is stationary when the fields p and p satisfy the Euler-Lagrange equations 

0(r - b 2 ) 

2p-2iy{p-p) = K 2 ' 

p 2 

and the condition of maximum for iA u must be determined independently for the various 
solutions by studying the stability of the action functional for arbitrary variations of the 
fields (cfr. sec.lgil). 

At the semiclassical level we are considering, it is assumed that the sum over field 
configurations represented by the path integral be dominated by just one solution among 
those satisfying the Euler-Lagrange equations fl38l) . Therefore, if p is such a solution 
yielding the dominant contribution to S, for symmetry reasons the solution p yielding 
the dominant contribution to must be equal to p*. In other words, we argue that only 
solutions with p(r) = p*(r) = pi + ip 2 are physically acceptable within the semiclassical 
approximation, pi j2 being real fields. 

Once the semiclassical solutions p = pi + ip 2 have been determined, the action is 
calculated by splitting the integration over r in eq. (157)1 into two pieces: (i) the first 
interval < r < b 2 where the fields freely evolve and the integrand is constant; (ii) the 
second interval b 2 < r < oo where the evolution is non-trivial but admits a constant of 
motion provided by the hamiltonian obtained by Legendre transform of the lagrangian in 
eq. J37J): 

H u = 2p lP2 + 2yp\ - -^— 2 = . (39) 

Pi + Pi 

No additional constant of motion exists, so that the system (138]) is not solvable analytically. 
Nevertheless, various relations [7J allow one to express the inclusive action in terms of few 
parameters of the solution, namely the slope t — t\ + it 2 of p(r) in the free region (i), and 
the asymptotic value p^ = p 2 (oo), according to the formula [7J 



__ 3 t 2 

Po ° 2t 2 + t 2 



(40) 



"2. 

This equation shows that the inclusive action vanishes for real solutions (p 2 = 0) in which 
case the 5-matrix is unitary. On the contrary, absorption is present only with a non- 
vanishing imaginary part p 2 . By studying the "inclusive" equation of motion (I3"8|) . we 
can determine for which values of b and y the solutions necessarily develop an imaginary 
part, and we can then compute the unitarity defect. However, since no analytic solution 
is available, we cannot exhibit an explicit criticality equation, and we have to resort to 
numerical or approximate methods. 



3.2 Real and complex inclusive solutions 

In the previous section we showed that a sufficient condition for unitarity is the existence of 
real solutions p = p* = p, because the inclusive action vanishes identically, thus implying 
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(0|S r ''S'|0) = 1 at semiclassical level. In this case both equations in (!38|) reduce to the 
equation ( II Op governing the elastic S-matrix, with the same boundary conditions, hence 
with real solutions only for b > b c . As a consequence, in this regime the model is unitary. 
This argument shows that the same critical point b c - - found to be b\ = ^^-R 2 for 
the elastic channel at y — — also governs the inelastic unitarity of the S'-matrix in a 
y-independent way. 

The above reasoning also implies that, below the critical impact parameter (b < b c ), 
only genuine complex solutions of (1551) exist. In order to characterize such solutions, we 
have written a numerical program that seeks for independent pairs (p, p) satisfying the 
inclusive system. We found that, without imposing the "physical" condition p = p*, there 
are four distinct complex solutions for any given value of b and y. By labelling each 
solution with the value of t = p(0), we can represent the four solutions with four points in 
the complex t-plane. At fixed y and increasing b these points move, spanning the curves 
Mi U M 3 , M 2 U M 4 , Ni U N 3 and N 2 U N 4 depicted in fig. [2k, in the direction of the arrows. 
Now, if we consider only the solutions with p = p*, their number varies with b according 
to the following scheme: 

• b < b c : there are two complex solutions, one with positive (Ni) and one with negative 
(N2) imaginary part; 

• b c < b < bd(y): there are two real (M3, M4) and two complex (iVi, N2) solutions, 
the latter with positive imaginary part; here 6<j is a ^/-dependent critical value that 
will be analytically described in the next section; d 

• b > bd(y)~- only the real solutions (M3, M4) survive. 




As already mentioned, we can accept only field configurations corresponding to max- 
ima of iA u . Analytically, stability is expressed by requiring the second variation of S^L 
to be positive definite with respect to arbitrary real variations (5p, 8p) of the solutionis 

5 At b — b c the three solutions Mi, M2 and N2 are pure real and coincide. At b = bd the two solutions 
Ni, and N 2 coincide. 

6 We consider only real variations of p and p for compatibility with the definition of the path inte- 
gral (I19[) which is defined by integration over real fields. 
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By definition, the first-order variation of the inclusive action vanishes on the solutions: 
5A U = 0. The imaginary part of the second-order variation yields 

5 2 3A = / [y(Sp - 5pf - 0(r - ^yj'^P [(M 2 + (Sp) 2 ] } dr , (41) 
which is positive definite provided 

p 2 (3p 2 1 -p 2 2 )<0. (42) 

In practice, when \p 2 \ is smaller than |pi|, as happens in the important region b ~ b c , 
stability occurs for p 2 < 0. Therefore, the only complex solution satisfying conditions f )42p 
is found on N 2 \^ t<0 f° r b < b c ; all other complex solutions are unstable. On the other 
hand, any real solution is stable. 

To summarize, for b > b c we have two real and y-independent acceptable solutions 
which coincide with the exclusive solutions studied in [6], providing a unitary ^-matrix. 
For b < b c instead, we have only one complex acceptable solution with negative imaginary 
part < 0) and corresponding negative action iA u < (cfr. fig. [2b), thus determining 

a unitarity defect for the S'-matrix. 



3.3 Behaviour around the critical point 



In order to investigate the onset of absorption below the critical point and to quantify 
its magnitude, we study the behaviour of the solutions of the inclusive equations (j3"81) for 
b ~ b c . In particular, we consider the physical (p = p*) complex solution N 2 which provides 
the unitarity defect for b < b c . The fact that this solution becomes real (p 2 {j) = 0) for 
b = b c , suggests to perform a perturbative analysis in which the imaginary part p 2 is 
considered a small quantity in some neighbourhood of b ~ b c . 

We first set up the relevant equations in a way that is convenient for the perturbative 
expansion. In terms of the real components pi ;2 , eqs. f)38p read 



2p 1 + 4yp 2 = Q(r-b^ Pl p2 



(Pi + Pl) 2 
2 -2 Pl p 2 



Pi(0) = 0, p 2 (0) = 0, Pl (oo) 



p 2 (ooj 



The evolution for r < b 2 is linear in r: 



PlA T ) = t l,2 T ) 



(r<b 2 



(43a) 

(43b) 
(43c) 

(44) 



For t > b 2 , thanks to the existence of the integral of motion fl39|) . we can reduce by one 
the order of the differential system fj4*3]) . e.g. replacing the system fj4*3|) with the linear 
combination fl45aD -2^ fl43b"l) and eq. (J3HD itself. 

A convenient way of rewriting eq. fl3T?|) is to consider pi the independent variable, 
i.e., p 2 = p 2 (pi(r)). This is possible since pi(r) turns out to be a monotonic increasing 
function of r. By denoting with a prime the derivative with respect to pi, we have 



P2 



dp2 
dpi 



P2 
Pl 



(45) 
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Dividing eq. (139]) by (pi) 2 and rearranging the factors, we end up with an equivalent form 
of the inclusive equations 



2pi 



Pl - p\ + %PlP2 

(Pl + Pl? 
1 



k P2 pi(pi + pI)(i + 2/p 2 ) 



(46a) 
(46b) 



which is particularly suited for the expansion in p 2 we are going to perform. 

The crucial observation is that, to each perturbative order in p 2 , there exists a second 
integral of motion which allows us to express p 2 as a function of p\ and of an "integration 
constant", e.g. p 2 (oo) = p^. This is easily seen by carrying out explicitly the calculation. 



Lowest order By expanding in p 2 the r.h.s. of eqs. (I46p at lowest relative order, we 
find for t > b 2 = b 2 



A- _i_ 

-2 2" 

p2 pi Pi 



(47a) 
(47b) 



Eq. (|47aj) is nothing but the exclusive equation f lTUj) which admits the integral of motion 



P\ + ~ = 1 

Pi 



yielding p\ as a function of pi, i.e. 



Pi 



1 . 

Pi 



(4f 



(49) 



Considering now eq. (147bj) . we can use eq. (T4"9"j) to replace pf in the r.h.s., so as to obtain 
the logarithmic derivative of p 2 with respect to pi and, after straightforward integration, 
p 2 itself: 



, dlogp 2 
dpi 



2^ 

P2 



2 2 
P2 = Poo 



1 

1 

Pl 



(50) 



where we have taken into account the fact that pi(oo) = oo and therefore the integration 
constant poo has to be identified with p 2 (oo). The matching at r = b 2 = b\ = 3v3/2 of 
the above solutions (I49f50l) with the free evolution (1441) provides the lowest order t- values 



9' 



(51) 



First order The first order relative corrections are obtained by expanding eqs. (146]) to 
first relative order in p 2 . 



2p 1 = l+4y£§ 
Pi Pi 



P2 



yp2 



pi pi 



(52a) 
(52b) 



13 



We now substitute the expression (!50|) into eq. (I52al) . obtaining 



1 1 / 1 d r . , 

2pi = — + Ayp^ — J 1 = — Vi(pi) , 

PI PiV Pi dpi 



whence the (first-order) conserved quantity 



(53) 



f,l = 1 - - V^OO)] = 1 - - + 1/Poo^ 

Pi 15 



If 13 

1-- 2 + 2 

Pi V Pi PI 



(54) 



where pi(oo) = 1 has been imposed. 

The first-order correction to p2 is now found by substituting eqs. (I50|54p into the r.h.s. 



of eq. (152bj) and, after integration in pi, we obtain 



pl = pI < i 



i 

pi 



+ ypc 



3pi 



1/11 28 

+ 



pi \l5p( 15pi 



(55) 



The values of t\ and t 2 are again found from the matching at r = b 2 with the free 
solution (I44p . However, at this level of accuracy, b is slightly different from b c and also 
t\ differs from t c . We parameterize these differences by introducing the adimensional 
parameters /3 and e such thatS 



1-/3 



U = Ul + e) 



(56) 



Note that /3 > means b > b c . Evaluating eqs. (I54ll39f55l) at r = 6 2 yields respectively 



/? = A iy t 2 , ^ = i(9 - 2^3) 



Poo — ~h 



I + — (27- Vs)yt 2 
15 



(57) 
(58) 
(59) 



To summarize the results so far, we have found the first-order corrections to t\ 
t c (l + e), t 2 and p^, obtaining 



f3 ~ e ~ yt 2 ~ ypoo • 



(60) 



/3 and e are of the same order, t 2 and p^ are of the same order but the latter have a 
y factor relative to the former. In order to deal also with the y = case, it is then 
convenient to write (3, e and p^ in terms of t 2 - 



Second order We can compute the second order quantities by expanding eqs. 



2pi 

l P2 



^5 



PI 



Pl 



,p3 
'pI 



1 / | 2/2 

.2 o v i - yp 2 + y P2 

Pl Pl 



pf 

p! 



(61a) 
(61b) 



7 This definition of (3 differs at O (/3 2 ) from the analogous definition of ref. [5] eq. (72). 
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Figure 3: [Left] The two solutions t 2 (red dashed and blue solid) versus (3 for y = 0.2. 
They coincide at b = bd- Only t 2 ^ for b < b c is physical (thick line). [Right] Comparison 
of the previous solutions for y = 0.2 (blue solid and dashed) with ^ in the y = 
case (green dotted and dash-dotted), the latter showing a square-root behaviour around the 
critical point (3 = 0. 



The explicit calculation is based on the strategy and the results presented at first order. 
Here we just write down the relations among the relevant parameters (3, p^, e and t 2 . 



(3 = A iy t 2 -^(l + A 2 y 2 ) t 



2 

2 i 



.4, 



291 - 112^ 
75 



1 , /-, if 16^3 - 23 „ i , 
--(12 - V3)yt 2 + - 1 + -1— y 2 ) t\ . 



9 

~t 2 
2 2 



2 , l 883 - 96^ „ , ., 
1 + —(27 - Vs)yt 2 +1-3 + I <2 



(62) 
(63) 

(64) 



We note the different behaviour of t 2 versus (3 at different values of y (fig. [3b). Without 
opened inelastic channels (y = 0) the imaginary component p 2 of the field grows in a non- 
analytic way as soon as b decreases below b c : □ 



~ T.n ~ 



'e . 



(65) 



On the other hand, with a finite inelastic emission (y > 0) the imaginary component 
p 2 grows linearly for small values of b c — b, i.e., \(3\ <C y 2 < 1, according to eq. flBTl) . 
In the intermediate region y 2 ~ <C 1, where both y and (3 are small, the quadratic 
equation f )62|) yields the two solutions (fig. [3^) 



94 ±} ^ Aij/ ± yAjy 2 - 18(3 (66) 

provided its discriminant is greater than zero, and this happens only if the impact pa- 
rameter b is smaller than a ^/-dependent critical value bd determined by 



A 2 2 

P < ^f- = PM , b 2 M = 



(67) 



3 Actually, the same square-root behaviour occurs in the region y 2 <C |/3| < 1. 
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The above condition reproduces the structure of the complex solutions numerically found 
in sec. 13. 2[ namely the existence of the two complex solutions iV\ and N 2 for b 2 < b 2 < b 2 d . 
For b > b c both have positive imaginary part (t 2 > 0)> while for b < b c only one of them 
(t 2 "*) becomes negative and is therefore physically acceptable according to the stability 
requirement. 

Such different behaviours affect the inclusive action. As a consequence, there are 
different trends of absorption, according to whether few graviton emission is allowed 
(y 0) or a finite contribution of the inelastic channels is considered. 



3.4 Unitarity defect 



The solutions just obtained with the perturbative method allow us to compute the uni- 
tarity defect e lAu — at least for small enough y and (3. In fact, by recalling eq. (PfOl) 



iA u = 4a 



Pc 



3 U 



2t( + t 2 



(68) 



the imaginary part of the unitarity action is given in terms of £1, t 2 and p^, whose 
dependence on the impact parameter b ~ b c has been analytically obtained in the previous 
section. 

Above the critical point b > b c , the physical inclusive solutions are real, P2(t) vanishes 
identically, and the same is true for t 2 , Poo and A u , thus implying a unitary S-matrix. 
On the contrary, below the critical point b < b c , the inclusive action is governed by the 
complex perturbative solution t 2 ~^ < which joins continuously with the real solution(s) 
at b = b c . In order to study the transition across the critical point, we expand the 
action fl68l) in t 2 by means of eqs. (163]) and fl64|) . where we recall that t\ — (1 + e)/v3. 

We note the interesting feature that the linear terms in t 2 cancel in the expansion of 
A u . Therefore, the suppression of the S'-matrix starts at second order in t 2 . 



i-A u 



a 



12 



(9 - 2V3)yt 2 2 + 18 



509 - 168^ 
75~ 



-y 



4 + 



(69) 



Due to the interplay between yt 2 and t%, A u is characterized by various regimes. 
At y = the action (I6"9"|) is cubic in t 2 ~ —^—[5: 



\A U 



a 



ist 3 2 ~ 



4V2 



3/2 



SAel 

a 



(70) 



and reproduces the fractional exponent 3/2 of the elastic action (fl6|) as shown in ref. [5]. 
Note that iA u is negative due to the choice of the stable solution t 2 ■* < 0. As a conse- 
quence, the S'-matrix is suppressed for b < b c . The same happens at small y \t 2 \ ~ 

At intermediate y ~ 1^1 <C 1, the two terms in eq. fIBTJl) are of the same order and the 
S- matrix is always suppressed for b < b c since t 2 < 0. 

At finite y ^> yfjsj, the term proportional to yt 2 dominates the inclusive action, which 
shows a negative-definite quadratic behaviour in b — b c , with a vanishing maximum at the 
critical point: 

a 5 A-i y 



(71) 
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This is in agreement with the numerical result represented by the solid line in fig. [2b. 

The above analysis for the unitarity action can be summarized by distinguishing three 
regimes, at least for y < 1: 

• For b < b c and (b — b c ) 2 ^> y, the unitarity action iA u is negative and decreases in 
size with a power-like behaviour (b c — 6) 3 / 2 as b — > b c . This includes the case y = 0. 

• For b < b c and (6 — b c ) 2 <C y, the unitarity action iA u is negative and vanishes 
quadratically at b = b c . 

• For b > b c the unitarity action vanishes identically. 

It is important to realize that the proper choice of the physical solution avoids a unitarity 
excess (0|S'^5'|0) > 1. Such an unphysical behaviour would occur if we don't reject the 
solution, since the corresponding inclusive action iA u (t^) would be positive for b < b c , 
as shown by the circles in fig. [2b. 




e-14 \L 



0.5 1 1.5 2 

b 

Figure 4: Numerical evaluation of the unitarity deficit (solid lines) at various values ofy, 
and comparison with the quantum v.e.v. squared of the S -matrix (dashed lines) illustrating 
the increasing contribution of the inelastic channels at larger y. 

The numerical evaluation of the ensuing S-matrix suppression for various values of y is 
presented in fig. SI where we plot the semiclassical estimate of the v.e.v. (0|S'^S'|0) = e 1 " 4 " 
and compare it with the n = vacuum contribution |(n = 1 -S' 1 0) | 2 of the unitarity 
sum (1341) . The results show that, if y is small, there is a considerable unitarity defect 
in the S-matrix, the inelastic channels providing only a small correction to the elastic 
suppression. At larger values of y there is an important recovery of unitarity due to 
the contribution of the inelastic channels. However, it must be kept in mind that y - 
related to the maximum rapidity of the emitted gravitons — cannot be arbitrarily large, 
because energy conservation limits the energy of the emitted particles to be smaller than 
the available energy y/s. Actually, energy conservation prevents y to assume large values, 
which at most can be of order y ~ O ( 1 ) . 

In conclusion, a unitary defect is always present in the semiclassical estimate of 
(0|S'^S'|0) when the impact parameter b is smaller than the critical one b c . In this re- 
gion the contribution of inelastic processes only partially compensates the suppression of 
the elastic channel, at least at semiclassical level. 
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4 Search for quantum transitions 



The semiclassical estimate of unitarity of the S-matrix performed in the previous sec- 
tion gave an exponential suppression of (0|ifv5'|0) ~ e~ Gs for b < b c . In this section we 
investigate whether quantum effects provide larger contributions eventually restoring uni- 
tarity. The quantity that we are going to study at quantum level is the eigenvalue of the 
S'-matrix operator (TH?|) . This is possible because the specific form of the S'-matrix — a 
coherent-state operator acting on the Fock space of gravitons. — allows us to determine 
its eigenvectors and eigenvalues [7], as we briefly recall. 

In order to construct such eigenvectors, it is convenient to introduce (normalized) 
graviton-coherent-states which are parameterized by arbitrary (complex) profile functions 
t](t) as follows: 

\ v (r)) = e -3fo*'">efo*' ot )|0) , (72) 

with the short-hand notation (rj, () = J °° i](t)((t) dr in terms of the azimuthally-averaged 
annihilation operators of eq. ([6]) 

a(r = * 2 ) = r^=^§ [a(r), a\r')] = 6(r - r') . (73) 

Jo ^V 71 " VY 

Since the action of the S-matrix operator on those states is a superposition of coherent 
states with shifted parameter 

S\ V ) = J [Dp] e~^ L ^ ^ \ v + i8 p ) , (5 P = v / 2^(l - pj) , (74) 
the functional Fourier transform of coherent states with imaginary parameter 

| Mr)}) = e i(^) / [V((t)} |iC(r)> (75) 



constitutes a complete set of S-matrix eigenstates. With the pre- factor e^' 1 ^ such states 
are normalized according to the delta functional 

(W(t)}\{uj(t)}) = 5({u - u'}) . (76) 

By acting on the states f T75|) with the S-matrix fTl9l) . it is straightforward to verify that 
the former are eigenstates of the latter. The corresponding eigenvalues, that we denote 
with e 1 " 4 ^, obviously depend on the real functional parameter w(r), and are given by the 
path-integral 

/■p(oo)=l 

eigenvJS) = e iA[oj] = / [Dp] e ^im+KS P ,^ _ ( 77 ) 

Jp{0)=0 

Endowed with eigenstates and eigenvalues of the S-matrix, we can perform a quanti- 
tative study of its unitarity properties at quantum level. For instance, we can reconsider 
the v.e.v. of S^S by inserting the completeness of u states, obtaining 

(OIS+SIO) = / [Vu\ (0|S t |{w})({u;}|S|0) = / [Vu\ e ~\(^) e -^M . ( 78 ) 



In the r.h.s. of the previous expression, we note two elements that determine the order of 
magnitude of the v.e.v.: the "density" of the intermediate state (w,o;) and the eigenvalue 
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(actually, its modulus squared). If we find \{co}) states such that (u,u>) = 0(1) and 
Qv4.[u;] is much smaller than the tunneling exponent, we have a case for an important 
quantum effect, reducing the unitarity defect. 

At semiclassical level, the eigenvalue ( 177]) is found from the contribution of p which 
makes the action 

A[u) = - 



Hp, P, t ) - \/2ay(l - p)w 



dr 



(79) 



stationary. The ensuing Euler-Lagrange equation is similar to that of the elastic ampli- 
tude (TlOi) . but with an additional term representing an "external force" proportional to 
do;/dr which depends on the eigenstate: 



P 



0( 



T 



2p 2 




w(r) 

a 



50) 



It is this external force that might help to cross the repulsive Coulomb barrier thus 
avoiding the exponential suppression for proper values of oj[j). 

Here, however, we want to estimate the eigenvalue ( 177)) at quantum level. In order to 
do that, we note that e lA ^ can be rewritten as a matrix element — in the Hilbert-space 
quantizing the ID system ( 180)) — of a proper evolution operator £4,(0, oo), in complete 



(81) 



analogy to the quantum analysis of the elastic amplitude of sec. 12.31 

e M M = (p = 0|Z4(0,oo)|p=l). 



The bra and ket states embody the boundary conditions of the functional integration, 
while the Moller operator governs the o;(r)-dependent dynamics, and is determined 
in the following way. Starting from the effective lagrangian defined by the integrand in 
eq. ( I79p . by Legendre transform we derive the conjugate momentum and the hamiltonian 



n 



dp 



a 



2a(p - 1) + y/2ayu(r) 
Q + 2a- y/Zayu\ 2 Q(r-b 2 ) 



2a 



P 



to 



(82) 



We quantize this system by imposing canonical commutation relations [p, II] = i, e.g., by 
identifying the operator II + 2a: = —id /dp in the p-coordinate representation. Finally, the 
evolution operator obeys the differential equation 



d 

i— Z4(t, oo) 



By defining the state 

\ip w (r)) ee Z4(r, oo)\p = 1) obeying 



(85) 



we can express the S-matrix eigenvalue ( 181)) as the amplitude of the wave function at 
p = and r = 0: 



*Pu(p,T) EE (/#u;(V)} 



Mb 



^(0,0). 



(86) 
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The above wave function can be determined by solving the Schrodinger equation (185b ) . 
Before doing that, we eliminate the w(r)-dependent shift in momentum and energy by 
the similarity transformation 

V> u (p, r) -> e 1 ^^ e iVT ^ dr ' # w (p, r) (87) 

so that the Schrodinger equation for \I> is 

1 <9 2 Gfr 



■ 9 T, \ 



*(p,r) 



The shift in momentum has produced a linear potential term due to a r-dependent external 
force, as already noted before eq. flHUj) . Since t/> and \I/ differ only by an unimportant phase 
factor, we can replace the former with the latter in eq. fl86|) and use the hamiltonian in 
eq. flHHl) to compute the time-evolution operator U^. 

Let us now compute \I/(0,0) in the case 6 = 0: in this way the r-dependence only 
comes from the external force u. We adopt the perturbative approach by splitting in 
an unperturbed term H u= o = H c which coincides with the "Coulomb" hamiltonian (j29l 
governing the elastic amplitude, and in a perturbation given by the time-dependent po- 
tential \f2ayujp. For the sake of simplicity, we assume that the external force is active 
only within a finite time range < r < T so that for r > T the evolution is just 
of Coulomb type. In this way the wave function can be determined by expanding the 
evolution operator in Dyson's series, yielding at first order 



* u (p, 0) = MP) - iVlay I dr w(r) (p\e iHcT pe 



|^c>H , (89) 



We obtain a more explicit expression by inserting a complete set of unperturbed eigen- 
states |0 (n) ) of energy E n before the operator p, by exploiting eq. ( |29|) . i.e., H c \ip c ) = 0, 
and then projecting on the position eigenstate (p = 0|: 

e u N| 6=0 ~ ^(0) - \^T y V^)(0) (<f>W\p\i>c) f T e iE ^u(T) dr . (90) 

We see that the effect of the interaction potential ~ up is to cause transitions from the 
initial state ip c towards other eigenstates of H c . Eq. (190]) can then be interpreted by 
saying that quantum effects provide a recovery of unitarity if the interaction potential 
induces transitions towards states with a wave function at the origin much larger than 
the tunneling amplitude: n (O) 3> V'c(O) ~ e~ na . We expect that such states could belong 
to two groups: 

• bound states whose wave function is significantly different from zero in the region 
p < and is finite at p = 0; 

• high-energy continuum states for which the suppression from the Coulomb barrier 
is smaller. 

More precisely, the entity of the contribution of those states to the eigenvalue e 1 " 4 ^ in 
eq. ( 190|) is given by three factors: the wave function at the origin n (O), the matrix element 
(4> n p\ipc) and the time integral. The first two will be studied in detail in the next sections. 
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As for the time integral, we need to specify the function u>(t). By analogy with the 
quantum mechanical phenomenon of induced transitions, we consider as a good candidate 
of external force an oscillatory function which can induce big transition amplitude at 
resonance. More precisely we choose 

u(r) = ^=sm(Kr)Q(T-T), (91) 



where the normalization has been chosen to keep (uj,u) = O (1) in order to avoid the 
exponential suppression in eq. fl75|) . Whith this choice, the time-factor reads 



' T K 

(92) 



T n = [ e iE - T u(r)dr 
Jo 



2VT 



i(E n -K)T _ i 

+ (K —K) 



L i(En-K) 



We are interested in situations where the perturbative corrections to the eigenvalue (1901 
are larger than the lowest order term. In this case, the square modulus of the eigenvalue 
can be written as 



| e u[< 



I^Vcc'TT^^Vc C* r^-E,,)T M(En ~ K)T\ sin[(E n , - K)T] 



(E n — K)(E n i — K)T 

(93) 

where the c n 's contain the p matrix elements and the wave function at the origin. For long 
interaction time T, the preceding expression is peaked for E n = K and provides a delta 
function of energy conservation K 2 5(E n — K) which suppresses the interference terms 
with respect to the squares of the individual amplitues. The latter provide a contribution 
of order 0(1) around the states compatible with conservation. Therefore, the time factor 
is 0(1) if (u,u) = 0(1). 



4.1 Transitions to bound states 

According to eq. (190]) . the entity of quantum effects to unitarity depends crucially on 
the matrix elements (<j) n \p\tp c ) and on the wave function at the origin 4> n (0) of the state 
reached by the induced transition. In this section we determine both factors for the case 
of transitions to bound states. 

In fact, the Coulomb potential ~ 1/p in the unperturbed hamiltonian H c is attractive 
for p < giving rise to a discrete spectrum associated to bound states which could 
correspond to collapsed states characterized by a strong gravitational field h oc (1 — p). 
The discrete spectrum of H c is given by the "energy" eigenvalues 



E n = -a - 



a 



in 



1)2 



[n 



e N) 



(94) 



and the (normalized) eigenfunctions are expressed in terms of the irregular Whittaker 
function [16] W (closely related to the confluent hypergeometric function U) 



MP) = N n 



T (l + n ) w i y i(x) e(x) +r(f-n) w n+u (-x) e(-x) 



2"2 



n+ 



2'2 



,(0) 



ca 



(n+i) 3 / 2 



x 



4a 2 p 



(95) 
(96) 



21 



where c ~ 0.45. These wave functions are significantly different from zero for p < 0, and 
in particular where the energy is larger than the potential (see fig. At variance with 
■0 C (O) ~ e _7ra , their value at the origin n (O) is not suppressed with a, but only by a power 
of n, as can be read from eq. fl96|) . For this reason, transitions to bound states might be 
important for the estimate of the eigenvalue (190]) . The crucial factor is then the matrix 
element 

/+oo 
<P* n (p)pMp)dp (97) 
-oo 

which can be estimated by using the WKB approximations for n and ip c . We use the 
notation n = n + | and divide the integration domain into 3 regions: 
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Figure 5: Plot of 'Ripe (dashed) and of 4> n (solid) for n = (thick) and n = 2 (thin), the 
last function being multiplied by a factor of 3. Here a = 1 is small in order not to have 
a huge suppression of the wave functions in the respective classically forbidden regions. 



(p < 0) For negative p, ip c oscillates with amplitude ~ e _7ra because of the tunneling 
suppression. The generic bound state eigenfunction oscillates in the interval — (n/a) 2 < 
p < 0, reaching its maximum value in the leftmost half-period, with amplitude ~ an -7 / 3 , 
and then decreases exponentially towards zero for p — > — oo. Therefore, the convergent 
oscillatory integrand is uniformly bounded by a~ 1 h~ 1 ^ 3 e~ na and the ensuing contribution 
to the matrix element is exponentially suppressed in a. 

(p > 1) In this region, ip c is an oscillating function with maximum amplitude ~ a 1//6 . 
On the other hand, \<f) n \ < iV„r(l + n)n- rn e- 2a2p/rn is strongly suppressed. It turns out 
that the product (p n ip c is oscillatory and exponentially suppressed in a 2 , and so is the 
contribution of this region to the matrix element ( 19T|) . 

(0 < p < 1) In the intermediate region \ip c \ increases from the value |V'c(0)| — e _7ra to 
values of order one, like e _7ra e 4a v / ^ ; while (j) n starts from (j) n (0) ~ a/n 3//2 and goes to zero 
as (i) N n e~ Aa ^ if p < (n/a) 2 ; (li) N n Y{\ + h)(x + n)-™ e -^ 2 p/™ if p > (h/ a ) 2 . It turns 
out that the product of wave functions 0„V> C is also in this case exponentially suppressed. 
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In fact, in case (i), the two esponential with ^/p cancel out and the product is of order 
e~ na ; in case (ii) we have 

\Mc\ < e-™N n VMe- h e f{x > h) , f(x, n) = ~ + 2v^ - n In (l + ?) , n = n + \ . 

(98) 

The exponent f(x, n) has a maximum at x = n, with value / max = (| — In 2)n. Therefore, 
the integrand is uniformly bounded by a (decreasing) function of n times ae~ na . Therefore 
also in this intermediate region the matrix element is exponentially suppressed in a. 

These results are confirmed by numerically computing the matrix element ( 19 7p . as 
shown in fig. El In conclusion, at least in this first-order perturbative treatment, the S- 
matrix eigenvalue ( 190]) does not receive significative enhancements from the bound-states 
eigenfunctions beyond the exponentially suppressed contribution ip c (0). 
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Figure 6: Numerical computation of the matrix element (9l\ ) for n up to 10. Each line 
represents a different value of a, here ranging from 2 to 6. 
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4.2 Transitions to continuum states 

The second mechanism of possible restoration of unitarity that we consider is the tran- 
sition towards "high-energy" eigenstates {4>e\ of the continuum spectrum of H c . The 
motivation is that, for p < 0, we expect the corresponding eigenfunctions <Pe{p) to be less 
suppressed than ip c , which is just the zero-energy eigenfunction. The eigenfunctions obey 
the Schrodinger equation {h = 1) 



H r <f>R = 



1 d 2 



4a dp 2 



+ a 



(99) 



and the condition of having only transmitted wave for p < 0. If we write the energy 
E = a (tg — 1) in terms of the parameter to > 0, it is easy to check that, with the 
rescalings a — > a/t and p — >■ pt^, eq. (1991 reduces to eq. (1291 for tp c = <f) . Therefore, the 
generic eigenfunction can be expressed in terms of ip c , namely 



>(to>;^) 



(100) 
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where the prefactor \fafn has been chosen in order to obtain the continuum normalization 

J <f>*M4>Eip)dp = 5(t -t' ). (101) 
Eq. fllOOp clearly shows that, for large energies E 3> a •<=>■ to 3> 1, the wave function 

net 

at the origin 0b (0) ~ e 'o is really much larger than ip c (0), because it is easier to cross 
the barrier at p > 0. 

It remains to evaluate the matrix elements (4>e\p\iP c )- By using the integral represen- 
tation of (fi c derived in ref. [7] and the relation f llOOp . we can write 



>E 



fa ,a, r +oc ft - e i2aptto 



77 fV 

e " 2 



N(a) = r- r— t: — r . (103) 

v ' cosh(vra)r(ia) v ' 

With this representation we can explicitly perform the p integration in the matrix element, 
as outlined in app. |A] The result is expressed in terms of hypergeometric functions: 

™ <t> E {p){4c{p) dp = N(a)N*(^)-^^I(t ) (104) 
jo to 4cra c^o 

ft — 1\ i<Q+Q °) r 

/(t ) = -5 — - e ™° vraoC coth(7ra)F(l - ia, 1 - ia ; 2; C) 

\t + 1/ L 

+ wT 1Q) v (1 ,~ 1 ?? F ( - ia > ~ ia ^ 1 - K« + ab); 1 - C)l , (105) 
r(l - i(a + a )j -I 

where C = -(T=^j2 and a = f . 

Let us now study the order of magnitude of /(to) which determines the matrix element 
to continuum states. We firstly consider transitions to high-energy states (£g = 1 + 
^ ^> 1). In this case, the variable £ ~ — ^ in /(to) tends to zero, so that the two 
hypergeometric functions assume finite values and cannot give exponentially enhanced 
contributions. Therefore, the order of magnitude of /(to) is determined by the factor e na ° 
in eq. (I105p . By taking into account the wave-function normalizations, the matrix element 
is of order 

N(a)N* f-)/(t ) T7 \ i / ^ ^ re™ ~ e"™ (106) 

to cosh(7ra) cosh(7rao)r(ia)r(iao) 

This result shows that the matrix element for transitions to continuum states of very high 
energy suffers a suppression comparable to that of tunneling. As a consequence, unitarity 
cannot be restored by these quantum effects. 

Secondly, we consider transitions to lower energy states characterized by to ~ 1 but 
a — ao — ^(to — 1) ^> 1, corresponding to a quite large jump in energy, implying a lower 
suppression of 0e(O) ~ e _7ra °. In this regime, £ — > — oo and the corresponding behaviour 
of the hypergeometric function is 

F(l - ia, 1 - ia ; 2; C) ~ e™ 
F( - ia, -ia ; 1 - i(a + « ); 1 - C) ~ A^^-^ >e~™ + A 2 e~™° . 
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The leading term is provided by the first hypergeometric function, yielding /(to) ~ e 2lTao . 
By inserting the wave-function normalization factors (of order e~ 7r ( a+a °' ) ), the tunneling 
amplitude 

M^eWc) ~ e -™°e-^ Q - Q °) ~ e~™ (107) 

turns out to be exponentially suppressed also in this case. 

In conclusion, the first-order perturbative contribution to the eigenvalue (190]) are expo- 
nentially suppressed by the matrix element in all those cases where we expected sizeable 
effects thanks to the enhancement of the eigenfuncions at the origin. Therefore, the 
quantum effects do not modify the semiclassical picture described in sec. El 



5 p(0)-fluctuations and short-distance singularities 

In the preceding sections we have definitely confirmed the unitarity defect previously 
found at semiclassical level for b <b c [7|: indeed, the critical behaviour of the action on 
the physical solutions around b ~ b c is as expected (sec. [3]) and the extra-contributions 
to the S- matrix elements due to quantum transitions are likewise suppressed (sec. H|). 
However, we have kept throughout the analysis the ACV boundary conditions p(oo) = 1 
and p(0) = 0. While the former is needed in order to match the large-distance behaviour 
with perturbative gravity and is thus unavoidable, the latter insures that the solutions 
being considered are UV-safe, so as to make the effective action self-sufficient. Is this 
really required? What happens if we let p(0) fluctuate? 

Let us recall that a non-vanishing value of p(0) = po (of order R 2 ) implies singularities 
at r 2 = of 0(r 2 ) and of h(r 2 ), as follows 

(2tt) 2 0~-^, h = V 2 ~ ~^(r 2 ) , (108) 

where the latter (implied by the outgoing flux of V0) can be interpreted by assigning 
p(r 2 ) a discontinuity p at r 2 = so as to have 

(irR) 2 h(r 2 ) = (1 - p) reg - Po R 2 6(r 2 ) . (109) 

Although the form of such singularities is probably model dependent, their existence is 
expected for generic boundary conditions. They stem from the very fact that R 2 is the 
dimensionful coupling of the 2-dimensional action and in fact they occurred already 
in the more general 2-dimensional treatment of [3]. 

The above singularities produce a singular metric (fT2"j) and a divergent action. If we 
take them seriously, the divergent part of the action in eq. ( f!9|) is 



poo 

A= -a 

Jo 



;i"P) 2 -ye(r 2 -6 2 ) 



r ] r 2 poo j 2 

— + V / 2^/ o (I ~ p)Hr 2 ) + a\r 2 )} — 

(110) 



= Aiiv + Ae g , Aiiv(po) - -aPo^(O) - ^2ayp S(0) (a + aj) > ( m ) 

where the formally divergent 5(0) ~ O is presumably regularized by the string in 

the ACV approach. Note that the divergence affects not only the c-number part of A, 
but also the r 2 = mode of the coherent state operator, that we have normalized to 
[ao,%] = 1- 
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At this point we notice that the divergent part of the S-matrix suppresses in a drastic 
way elastic and quasi-elastic processes with p Q ^ because 

2 

(0| e u div(po)|o) ~ e -a(i+2/)p^(o) ^ e -«a+w)^% ( 112 ) 

shows violent oscillations and absorption, with exponent of order ~ ctp^jr ~ ct 2 Po f° r 

A s ~ Ap = V G/i. Therefore, at quasi-elastic level, we were justified in setting p = in 
the first place. And that would close the argument. 

The above conclusion is not fully satisfactory, though. On one hand, p Q ^ is associ- 
ated to a singular metric that we think corresponds to classically trapped solutions [T7H2T], 
and we would like to know about their fate at quantum level. On the other hand, if we 
think of the action (lllOp from the unitarity point of view, the divergent part is, after all, 
a hermitian operator and we expect it to contribute a unitary S'-matrix if we sum over all 
possible emissions. This statement is confirmed by the inclusive equations (sec. [3]) because 
we do have real- valued solutions for b < b c also, provided we take p = p = p m ip) > 0, 
as mentioned at the end of sec. I2.2L In this cases, the inclusive action vanishes: should 
we conclude that considering po ^ insures S'-matrix unitarity of the model? 

The trouble with our second argument is that it is inconsistent with energy conserva- 
tion. In fact, the action (11101) provides total probabilities of order unity only in association 
with quite a number ~ a-^- of hard emitted gravitons with energies of the order of the 
Planck mass. Barring other dynamical effects (possibly providing a red-shift), this would 
require an energy which is much larger than yfs. 

Nevertheless, from the above considerations, our tentative conclusion is that, from 
the unitarity stand point, we should actually consider the solutions with p 7^ also by 
keeping in mind that the present model is inadequate in the short-distance region and 
we should therefore come back to the string dynamics, that was originally neglected in 
the regime R 3> X s . The role of the latter is not only that of regularizing the would- 
be singularities at the string scale, but also of providing a number of physical effects at 
distances which are intermediate between X s and R. 

For instance, a well-known [UI22] effect, occurring already at distances of order R, 
is string excitation by tidal forces (or "diffractive" string excitation). It looks straight- 
forward to estimate it in the present model, but it has not been done yet. Another 
issue, perhaps more important for the unitarity problem, is the evaluation of rescatter- 
ing corrections [I]. Here the string is needed in order to regulate the UV divergencies 
in the longitudinal variables so as to extract a finite answer. We expect from them an 
improvement of the longitudinal dynamic^] for b = O (R) and perhaps some hint as to 
the existence of trapped solutions for p ^ 0. If effects of this kind turn out to yield 
an important contribution to unitarity, then this will involve relatively soft gravitons, of 
energies ~ h/R and, for this reason, may be consistent with energy conservation. 

It may be also that string effects are mostly confined at scale X s . In such cases we 
expect that compactified dimensions and/or universes can be excited and, therefore, that 
probability — in our original scattering process — is lost for good reasons: in order to 
recover it we would need scattering data in some extra worldo Whichever the case, we 

9 For instance, we expect the wave-fronts of the solution (TTJ) to be shifted [23J when they merge each 
other just to make their motion consistent with the scattering process. 

10 An example of this kind is perhaps provided by the recent analysis of string-brane scattering [24] in 
the case that the metric potential is singular enough to allow fall into the center. 
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are inclined to think that, for b < b c , solutions with p(0) 7^ 0, depending on the string 
dynamics, can play an important role in explaining the unitarity loss of the present model. 

6 Discussion 

The first and most detailed outcome of this paper is that the reduced-action model for 
the transplanckian S- matrix, valid in the regime R,b ^> A s , shows a unitarity suppression 
in the region b < b c , that we think corresponds to classical gravitational collapse. This 
feature, originally found at semiclassical level in [7], is here confirmed at more general 
quantum level in sec. H]and by a perturbative method around the critical region in sec. |3] 
We find in particular that the unitarity defect around the critical point is non analytic 
with a fractional critical exponent when inelastic emission is forbidden (y = 0), while it 
shows a more regular behaviour when y > 0. 

Since such unitarity defect raises questions about the information loss in gravitational 
collapse, we have tried to investigate also whether the present model is complete and if 
not, how to complete it. Let us first note that the model has a hermitian lagrangian and is 
thus expected to be unitary, barring some special conditions. But a key role is played by 
the boundary condition p(0) = 0. On one hand, it is needed in order to provide UV-safe 
solutions with regular metric, so as to insure that the eikonal is a perturbative series in 
R 2 /b 2 for b > b c , independently of the string scale A s . On the other hand, that condition 
is the main cause of the unitarity defect for b < b c because, in order to reach p = at 
r 2 = 0, the amplitude has to cross a potential barrier in p-space by a tunnel effect, and is 
thereby suppressed. 

For the above reasons, here we suggest to consider the solutions with p(0) 7^ as the 
alternative path that the probability flow may take. Such solutions yield a singular metric 
and a divergent action and probably correspond to classically trapped fields for b < b c . In 
the toy model of sec. |S]we let p(0) fluctuate with a weight provided by the reduced action 
itself, and we find that p(0) 7^ is violently suppressed at elastic and quasi-elastic level, 
thus justifying the condition p(0) = in those cases. Nevertheless, the model is formally 
unitary and acquires a total probability of order unity in association with a large number 
of hard gravitons, a process which is, however, forbidden by energy conservation. 

The paradoxical features above illustrate the problem we have in looking at the p(0) 7^ 
solutions: the present model suggests that unitarity might be recovered at inelastic 
level, but is by itself inadequate. Therefore, we have to rely on the string dynamics and 
related vertices pJJ] in order to describe the short-distance behaviour of the solutions 
and their contributions to unitarity. A key point is to determine the important scales of 
the evolution from R to X s , and the kind of inelastic channels which are involved, which 
may include states propagating in extra dimensions or universes. In sec. E]we have listed 
several string effects, some of which at scale X s and some at scales of order R, but only a 
careful analysis can tell us where the lost probability goes. 
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A Matrix element to continuum states 



In this appendix we explicitly perform the computation of the matrix element to con- 
tinuum states of eq. (I104p . Starting from the integral representation (11021) for (pE and 
ip c = 0o we have 



/ (j) E p^ c dp = N(a)N*(-) / dtdpdp 

J -co '0 J— oo 



x pe i2 ^— ) (t^lY ( ^Y t0 Siga }j ~ ~ 1} . (113) 



t+lj \P+1J (t 2 - l)(p 2 - 1) 

The integration in p can be performed by introducing cto = f- and rewriting the integrand 

as 

1 a J2a P (t~t 0P ) ( t~l \ " ( P^±\ ~ 1Q ° sign(t - l)sign(p - 1) 
2ap 0t ° e \t + l) \p + l) (t 2 -l)(p 2 -l) ' 1 } 

It produces the delta function S[2a(t — top)) which is used to perform the integration in 
t. Apart from the N normalization factors, the matrix element is now 

m_ r f t Q p - 1 V" ( p - 1 \ ~ 1Q ° sign(t p - l)sign(p - 1) dp 
2a 2 to J \t p+lj \p+lj ^-1)^-1) p - K ) 

By rescaling the variable top — > p the four branch points are found at p — ±£ , ±1 and 
the integrand gets a factor i 2 , that we conveniently rewrite as [£q — p 2 } + p 2 , obtaining 

ia 



/rpr-FT iTT „ /" /p - 1 Y" (v - to \ ™° sign(p - l)sign(p - t ) dp 



ivr f (p - 1 \ f p-t \ sign(p - l)sign(p - t 



+ 2^ s *'A?nJ li^J (p>-d(p»-« pdp - (116) 

Let us now consider the two integrals separately. If we perform the to derivative in the 
first integral, we obtain 

ivr f (p- \\ a fp-t \~ ia ° sign(p- l)sign(p - t ) dp 



2a 2 J \p+lj \p + t J {p 2 - l){p 2 - t 2 ) 

which vanishes, being proportional to the scalar product (4>e\4>q) °f two eigenfunctions 
with different energy. Only the second integral remains, and it can be simplified by noting 
that the factor p/(p 2 — tfy can be obtained by a further derivative with respect to to- In 
practice, 



4a 2 ao 

V ~ 1 \ ia ( P~t \ ~ ia ° sign(p - l )sign(p - t ) 
p+l) \p + t ) p 2 



I{to)= I (^=4) (f=^) ~ - " 0) dp . (118) 



The integrand of /(to) is built by four powers whose exponents' sum is —2. This allows 
us to express it in terms of the hypergeometric function, as follows. By a linear change 
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of integration variable x = ^Y a ^z^, we can map three of the four branch points into the 
"canonical" ones 0, 1, oo: 



i(a+oo) r+oo 



sign fx — 11 



x — 1 



x 



iao 



dx 



x(l — x) 



(119) 



where £ = ^_ 4 / o ° 2 ■ Here the integration path is slightly shifted off the real axis so as to lie 
below the cut [0, 1] and above the cut from — 1/C to infinity. The sign in the integrand 
amounts to split the integration into two pieces: one from — oo to 1 and one from +oo to 
1. By rotating the latter around 1 in counterclockwise direction, both pieces ranges from 
— oo to 1, one below and one above the [0, 1] cut. In this way, we arrive at the expression 



1 (t - 1 



2 V t + 1 



i(a+ao) 



7rao 



-2 cosh(7ra) 



x- iQ - x (l - xf^il - (x) iao dx 



+2 



x 



— icr— 1 / 



(1 + x) 1Q - i (l + Crr° dx 



(120) 



where the 2 cosh(7ra;) factor arises from evaluating the integrand above and below the 
[0, 1] cut. With the substitution y = in the second integral, we recognize in eq. (I120p 
two integral representations of the hypergeometric function, yielding 



Hto) 



to -I 
to + 1 



i(«+«o) 



7TQ0 



noioC coth(7ra)F(l — ia, 1 — i«o; 2; C) 



r(— ia)r(l — ian) / >i , . 

+ n •( " ia ' ~ iao; 1 - K« + «o); 1 - C ( 121 ) 

r(l - i(a + a )J J 



as in eq. f llOSp . 
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